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SUMMARY 

A method of analysis is developed- for two-dimensional flow on . 
general surfaces of revolution in turbomachines with arbitrary blade 
shapes. The method of analysis is developed for steady* compressible, 
nonviscous, irrotational flow that is assumed uniform normal to the 
surfaces of revolution. Incompressible solutions on a mean surface of 
revolution between the. hub and shroud are presented for four flow rates 
through each of two centrifugal impellers with the same hub- shroud con- 
tours but with different blade spacings. In addition, correlation 
equations are developed whereby the velocity components and the stream 
function distribution can be predicted for compressible or incompress- 
ible flow in straight-bladed impellers only, with any tip speed, flow 
rate, area variation, blade spacing, and for any flow surface of revo- 
lution. 


INTRODUCTION 

In order to achieve substantial improvements in compressor and tur- 
bine efficiency, new methods of design must be developed, and these 
methods should be based on detailed knowledge of flow conditions within 
turbomachines (compressors and turbines). The purpose of the analysis 
method, developed at the NACA Lewis laboratory .and presented in this 
report is to provide a tool whereby increased knowledge of the flow 
characteristics in turbomachines can be acquired. In addition the two- 
dimensional numerical results presented in this report provide the basis 
for future comparison .with the results of three-dimensional solutions. 

Several two-dimensional methods of analysis have r already beep 
developed (references! and 2, for example).. In general, these, varipup 
methods differ in technique and scope of application. In particular, 
in reference 2 a method of analysis., is developed for two-dimensional 
compressible flow in turbomachines with , conic flow surf aGes generated . 
about the axis of the turbomachine by straight lines. In this report , 
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the analysis method of reference 2 is extended to include arbitrary sur- 
faces of revolution generated about axes of turbomachines. In addition 
correlation equations are developed for the rapid estimation of flow 
conditions in im pellers with straight blades. 


MEfflOB OF AHALTSIS 
Ereliminary Considerations 

For given compressor operating conditions, a complete solution for 
the flow through turbomachines depends on the machine geometry (three- 
dimensional effects) and on the fluid properties (compressibility and 
viscosity) . The effects of compressibility are considered in the analysis 
method of this report but the effects of viscosity are not. Same three- 
dimensional effects are indicated by the combination of two-dimensional 
solutions. The idea of these quasi-three-dimensional solutions is con- 
sidered next. 

Quasi- three-dimensional solutions. - Consider the flow of an ideal 
fluid through a typical passage between blades such as shown in figure 1. 
The fluid is free to follow whatever path the pressure and inertia forces 
require of it. If, however, the -number of blades in the turbomachine is 
considered to approach infinity, the space between blades, as well as the 
blade thickness, approaches zero and the path of the fluid is limited to 
the curved mean surface of the blade. The fluid motion is thus reduced 
from a general three-dimensional motion to a two-dimensional motion 
restricted to the mean blade surface. The streamlines of this two- 
dimensional motion can be projected on the meridional plane (axial-radial 
plane) as shown in figure 2. The solutions are axially symmetrical 
because the infinite number of blades prevent the fluid properties from 
varying about the axis of the turb (machine. 

Axial- symmetry solutions establish mean circumferential flow con- 
ditions, but for finite blade spacing flow conditions vary circumferen- 
tially, as well as from hub to shroud. In order to investigate the cir- 
cumferential variation, which is superimposed on the circumferential 
mean flow conditions obtained by the axial- symmetry solution, it is 
assumed that the flow is limited to surfaces of revolution generated by 
rotating the meridional streamlines of the axial- symmetry solution 
(fig. 2) about the axis of the turbomachine. If these meridional stream- 
lines are selected with sufficiently close spacing, flow conditions can 
be considered uniform normal to these surfaces of revolution. Thus the 
flow is limited to two-dimensional motion on flow surfaces of revolution. 

Blade- to-blade solutions can be obtained for every flow surface of 
revolution generated by the streamlines in the meridional plane. There- 
fore, flow conditions can be determined approximately throughout the 
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passage. This resulting quasi- three-dimensional solution is obtained by 
the combination of two types of two-dimensional solution: axial- symmetry 
solution and blade-to-blade solution. Although this combination of two- 
dimensional solutions does not result in a rigorous three-dimensional 
solution, because, for example, the flow has been constrained to surfaces 
of revolution, the combination must result in a better picture of the 
true flow conditions than does either of the two-dimensional solutions 
alone. Analysis methods for the axial- symmetry type of solution are 
developed in references 3 and 4, for example. In this report analysis 
methods are developed far the blade-to-blade type of solution on flow 
surfaces of revolution generated around the axis of the turbomachine by 
known, or assumed, streamlines in the meridional plane. 

Coordinates . - The cylindrical coordinates R, 0, and Z shown 
in figure 3 are convenient to -use in theoretical studies of the flow in 
turbomachines. (All symbols are defined in the appendix.) These coor- 
dinates are dimensionless, the li n ear coordinates R and Z having 
been divided by the blade-tip radius rrp (so that -R equals 1.0 at the 

blade tip, for example) . The coordinate system rotates with the blade 
row about the Z-axis. The angular velocity a> is always considered 
positive and in the counterclockwise direction as shown in figure 3. 

Velocity components. - The velocity Q relative to the rotating 
coordinate system has components C^, and Q z in the R, 0, and 

Z directions, respectively (fig. 3). These velocities are dimensionless, 
having been divided by the stagnation speed of sound c Q upstream of 

the blade row, where 

c 0 2 = rgR*T 0 (i) 

In equation (1) R* is the gas constant, y is the ratio of specific 
heats, and T Q is the stagnation (total) temperature upstream of the 

blade row. The blade-tip speed is likewise dimensionless and becomes the 
conventional blade-tip Mach number Mp, which is defined by 

Mp = — (2) 

o 

Thus the tangential velocity of the blade at any radius R is equal to 
RMj nnd the absolute tangential velocity of the fluid is equal to 

(RHp + Qg) . ' 
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In addition to the three components 
Qg, and it is convenient to define 

lying in the meridional plane (RZ-plane) 
velocity Q M , where from figure 3 

2 2 

% = + 

The meridional velocity Qjj helps define the angle a, shown in fig- 
tire 3, from which figure 

Or = sin a (4) 

Qg = Qm cOS a, (5) 

Thermodynamic relations. - From the general energy equation, in the 
absence of heat transfer, it can he shown that along a streamline the 
static (stream) tengperature T is related to the relative velocity Q 
hy (reference 2) 


of the relative velocity Q^, 
a fourth velocity component, 
and defined as the meridional 

Qz 2 (3) 




(B %) 2 - Q 2 



( 6 ) 


where the subscript U refers to conditions far upstream of the blade 
row and where A is the whirl ratio defined as the absolute moment of 
momentum divided by r^c^ and given by 


X = R(HMu + Qg) (7) 

For isentropic flow the pressure ratio P and the density ratio p/p Q 
are likewise related to the relative velocity Q by 

r r 

p = - (^) T - j 1 + *jr [(E»t ) 2 - « 2 - T (s) 

and 

£“(i ) Tl= { 1 + I r [(k»d ) 2 - a 2 - T 1 (9) 

where p snfl p are the static pressure and density, respectively. 
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Analysis 

A streamline that was obtained from, an axial- symmetry solution in 
the meridional plane (fig. 2, for example) is shown in figure 4 together 
with the coordinates, velocities, and so forth at a point on the stream- 
line. At this point a tangent to the streamline in the meridional plane 
intersects the axis of the turbomachine with the angle a, as indicated. 
If this tangent line is rotated about the axis, it generates the surface 
of a cone that is tangent to the flow surface of revolution obtained by 
rotating the streamline (fig, 4). At the circumferential line of 
tangency a fluid particle on the flow surface of revolution is also on 
the conic flow surface. This fluid particle is shown in figure 5 on a 
developed view of the cone surf ape. The height ratio H of the fluid 
particle is defined by 

B - s; < 10 > 

where h (fig. 2) is the local spacing between center lines of the chan- 
nels formed by the adjacent^ streamlines (fig, 2) on either side of the 
streamline that generates the surface of revolution. 

Absolute irrotational motion. ~ In the absence of shock, viscous 
effects, and so forth, the absolute motion of the fluid everywhere on the 
surface of revolution is isentropic and hence irrotational so that the 
absolute circulation dT around every fluid partical is zero and from 
figure 5 


0T - 0 



Q?)B d6 


.1 a 

Qm«1 

y ® - w 

d 1 

to 

t 


da 


from which 


^0 

0 » 2R&J, sin a, + sin a + R sin a 

Continuity. - From continuity considerations for 
in figure 5 




(u) 


the fluid particle 


a 

m 



% 


_5_) 

sin «/ 


0 


( 12 ) 


where H is a function of R along the surface of revolution 

S - R(B) 
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A stream function if satisfies the continuity equation (12) if defined 
as 




If “£ Qm 


HR 


« - Qe 

Po ^ 


H 


sin a 


(13) 




Transformation of coordinates. - Eor convenience in the relaxation 
solution the following transformation of coordinates 1 b introduced: 


and 




1 dR 

(A0) sin"a R 


(14a) 


dij = 


d£ 

AO 


(14b) 


where AB is the angular spacing of the blades about the axis of the 
turbomachine. Equation (14a) can be integrated to give £ as a function 
of R with | arbitrarily equal to zero when R equals one, and equa- 

6 

tion (14b) can be integrated to give n equal to — ., By this trans- 
formation, which is conformal, the curved flow surface of revolution 
becomes a flat plane (Jq-plane) . Lines of constant £ and ti are shown 
transformed onto a flow surface of revolution in figure 6. The problem 
is solved in the |T]-plane and transformed back to the curved flow surface 
of revolution. The advantages of this transformation are: first, a 

warping of the blade shape which for straight blades results in straight 
parallel boundaries (blades) in the transformed plane, thus permitting 
equal grid spacing for solution by relaxation methods; and second, a 
reduction in the number of terms in the final differential equation for 
the stream function distribution, thus reducing the labor involved in 
obtaining a relaxation solution. 


G eneral differentia l e quation for stream function distrib ution.! - ” 
In terms of the transformed coordinates, equation (11) for absolute irro- 
tational motion becomes 


0 = 2RMji AB sin a + Qg A0 sin a + 
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and, equation (13) for the stream, function becomes 


= ^ Q M m (£9) 

If = “ ^ 


(16) 


Equations (15) and (16) combine to give (note that H »nd R are inde- 
pendent of 1}) 


p o 3. In 3 In ■£— 

5 ijr 3 ijr ^ In H 3ijr Po 3\Ir Po 3\Ir 

+ ^2 " H “ ^t M ~ Sn it| “ 


2 -fi. HR 2 (A0) 2 lip sin a = 0 (17) 

Po 

p 

In this equation the density ratio ~ is related to the velocity com- 

Po 

2 2 2 

ponents ^ and by equation (9) where Q = Q M + Qg , and the 

velocity components are related to the stream function ijr by equa- 
tion (16). This system of equations ((9), (16) , and (17)) can be solved 
for given boundary conditions by relaxation methods (references 2 and 5) 
to give the distribution of ijr for two-dimensional flow past arbitrary 
blade ‘shapes on . the flow surface of revolution. The velocity components, 
and therefore the pressure, density, and so forth, can be obtained from 
the distribution of ijr by means of equation (16) . 


Incompressible flow. - For incompressible flow the stagnation speed 
of sound contained in the definition of Q and % is undefined and 

■ Q 

its use can be eliminated by considering the ratio which is the 

Dam 

ratio of the relative velocity to the impeller tip speed. Thus, with 

-£• equal to one for incompressible flow equation (16) becomes 
Po 


d*i % 

ST “ % m ( 1 


_ % 


HR (A 0) 


(16a) 
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■where from equations (16) ^ is related to by 


*i = 


JL 

% 


so that equation (17) becomes after dividing "by Mrp 


d | 2 + a n 2 " 


2BR 2 (A0) 2 sin a = 0 


(17a) 


which is the partial differential equation used for the numerical solu- 
tions presented in this report. 


Method of Solution 

Numerical methods. - Solutions of the differential equations (17) 
or (17 a) for given blade shapes on given, arbitrary surfaces of revolu- 
tion and for given flow rates and "blade-tip ftfech numbers are obtained by 
relaxation methods using techniques discussed in reference 2. 

Boundary considerations. - As in reference 2 it can be shown from 
continuity considerations that if ^ is zero along the driving face of 
the blade (blade surface in the direction of rotation) the stream func- 
tion ijr-j- along the trailing face is constant and equal to 

= <P(A0) 

In equation (IB) the channel flow coefficient cp is defined by 


^o a T c o 

where W is the flow rate through the channel and is the channel 

flow area at the blade tip 

arp - r^h^Ae) (19a) 



.(18) 


For incompressible flow 
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The Joukowski condition and other "boundary conditions are discussed 
in reference 2. 

HUMERICAL EXAMPLES 

(NJ 

Solutions are presented for four flow rates through each of two 
centrifugal impellers with the same hub-shroud contours but with differ- 
ent blade spacings-. The solutions are for incompressible flow on a sur- 
face of resolution generated by the mean streamline of an axial- symmetry 
solution. 


Design and Operating Conditions 

Impeller geometry. - The hub- shroud profile of the impellers is 
described in figure 7. The impeller blades are straight and thin, and 
for mathematical simplicity extend indefinitely far upstream parallel 
with the axis of the Impeller. The diffuser is vaneless. Solutions 
were obtained for two blade spacings that, after accounting for the 
blade thickness, measured 14.67° and 32.80° about the axis. 

Operating conditions. - The fluid was considered to be incompress- 
ible and nonviscous. The solutions were obtained on a flow surface of 
revolution generated by the mean (0.5) streamline of the axial-symmetry 
solution in figure 2 (the R,Z- coordinates of which streamline were 
obtained from results of reference 4 and are given in table I together 
with sin a. and H) . Four solutions for different flow rates were 
obtained for each of the two blade spacings. At low flow rates an eddy 
(reference 2) exists on the driving face of the blade and at zero flow 
(closed throttle but impeller rotating) the eddy occupies the entire pas- 
sage. As indicated in reference 2 the eddy flow, which consists of fluid 
that remains in the impeller and rotates relative to the impeller with an 
angular velocity equal and opposite to that of the impeller, results 
from the condition of absolute irrotational fluid motion. The four 
solutions for each blade spacing were obtained for the following flow 
conditions: (l) eddy flow (no through flow so that the eddy occupies 

entire channel) , (2) zero eddy flow (through flow rate just sufficient 
to eliminate eddy) , (3) half eddy flow (half of the through flow rate 
required to eliminate eddy), and (4) through flow rate used in axial- 
symmetry solution (fig. 2 ) of reference 4. These flow rates can be 
expressed by ratios of upstream axial velocity to impeller tip speed 
(Q^/Mj^y. The through flow rate required to eliminate the eddy varies 

with the blade spacing so that the values of - for the four 

flow conditions described are given for both blade spacings in the 
following table: 



10 


HACA IN 2654 



f—) 

\%/U 

^''\Made spacing 
Flow^'"'~'^A0, deg 
condition 

14.67 

32.80 

Eddy flow 

0 

0 

One half eddy flow 

.060 

.105 

Zero eddy flow 

.120 

.210 

Flow rate of 

.3429 

.3429 

axial symmetry 
solution 




The solutions far the various ratios of flow rate to Imp eller tip 
speed were obtained by superposition of the eddy flow solution (rotating 
impeller with no through flow) with various percentages of the through 
flow solution (flow through nonrotating impeller) . 


Results 

The results of the numerical examples are presented in figures 8 
to U. The results are plotted on the surface of revolution between two 
impeller blades. The coordinate system is shown on the surface of 

revolution. The relation between |(AS) and the specified cylindrical 
coordinates R and Z along the surface of revolution was obtained 
from a numerical integration of equation (14a) and is given in, table I. 
From equation (14b) q is related to 0 by 



Streamlines , - The streamline patterns on the surface of revolution 
are shown in figures 8 and 9 for the 32.80° and the 14.67° angular blade 
spacingBj respectively. For each blade spacing the streamlines are given 
for the four specified flow rates. The impeller rotates in the clockwise 
direction so that the driving face of a blade is on the left of the 
channel and the trailing face on the right. Except for the complete eddy 
flow solutions (figs. 8(a) and 9(a)) the streamlines are designated by a 
stream function ratio such that the value of this ratio indi- 

cates the percentage of total channel flow between the streamline and the 
driving face of the blade. (For the complete eddy flow solutions the 
stream function ratio ^/(ijr^)^^ indicates the percentage of total 

flow rate circulating within the channel between the streamline and the 
nearest blade.) At each radius the spacing of the streamline is indi- 
cative of the relative velocity, with close spacing indicating high 
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velocities and wide spacing indicating low velocities. Thus, the veloc- 
ities are higher on the trailing face of the "blade than on the driving 
face, except at the "blade tip where the velocities are equal on "both 
faces and except for the complete eddy flow solutions (figs. 8(a) 
and 9(a)) where the velocities are equal on "both faces at equal radii. 

As indicated in the figures the fluid within the eddy remains in 
the impeller and rotates in the opposite direction to that of the Imp eller 
so that the velocities are inward along the driving face. (Also see 
reference 2, for example.) The eddy flow, which results from the condi- 
tion of absolute irrotational motion, is associated with changes in 
radius along the rotating flow surface and disappears as the upstream 
surface becomes cylindrical (constant radius). This situation is most 
clearly evident in the solutions for no through flow (figs. 8(a) 
and 9(a)). Thus, for quasi- three-dimensional solutions the eddy flow 
disappears when the surfaces of revolution became cylindrical. (This 
condition is commonly assumed when two-dimensional cascade data are 
applied to the design of axial-flow turbomachines.) For rigorous three- 
dimensional solutions the eddy continues to exist but its plane is normal 
to the axis of the turbomachine and results only in a distortion of the 1 
initially cylindrical surface. If this distortion is neglected by 
assuming the flow is constrained to surfaces of revolution, the eddy 
disappears on cylindrical surfaces. 

In the solutions for no through flow (figs. 8(a) and 9(a)) the 
greater magnitude of (^) miT1 for A9 equal to 32.80° than for. A 9 

equal to 14.67° indicates that the magnitude of the eddy flow rate 
increases with increasing blade spacing. 


Relative velocity components ■ - Lines of constant relative velocity 
.components are shown in figures 10 and 11 for the 32.80° and 14.67° angu- 
lar blade spacings, respectively. The meridional and relative tangential 

are plotted and only the solutions 


velocity components 


and 


% 


) 


for the eddy flow are considered. (For solutions with through flow. 


Q, 


9 


— remains the same and. 

Mrp 

that is proportional to 


% 


(©, 


Is increased an amount, given in table I, 
and is a function of £(A0) only.) 


% 

Lines of constant relative tangential velocity ^ are plotted in 

Q Q 

figures 10(a) and. 11(a) . The values of near the impeller tip are 
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negative so that the slip factor p (defined as the ratio of the average 
absolute tangential velocity at the impeller tip to the impeller tip 
speed) is less than 1.0. The resulting slip factors for the two blade 
spacings are 


Blade 

•Slip 

spacing 

factor 


V- 

(deg) 


14.67 

0.9111 

32.80 

.8142 


Qq 

Because — is unaffected by the through flow the values of p are 
■Mm 

independent of the flow rate (also see reference 6) . The values of 

Qe 


% 


at lower values of R upstream of the impeller tip approach a max- 


imum positive value and then decrease to zero as the surface of revolu- 
tion becomes cylindrical. It is noted that the maximum positive value 

Q e 

of — is considerably larger for the wider blade spacing. 


Si 

Lines of constant meridional velocity |j- are plotted in f ig- 

Qw 

ures 10(b) and U(b) . The values of — are negative along the driving 

±\*ti 

face and positive along the trailing face with zero values along the 

Qm 

center line of the channel between blades . The maximum values of -r— 

Mrp 

occur on the blade surfaces near the impeller tip and as the radius 
decreases -rr~ decreases and approaches zero as the upstream surface of 


Q, 


■e 


revolution becomes cylindrical. It is noted, as for tj-, 
blade spacing increases the magnitude of tt- increases. 


that as the 



2466 


HACA TN 2654 


13 


CQRREIATIOU OF RESULTS 

The relaxation solutions presented in this report required consid- 
erable labor and it is therefore desirable to have a more rapid method 
of analysis to serve as a useful engineering tool. The correlation 
equations to be presented in this section have been developed for this 
objective. These equations are limited to impellers with straight blades 
on arbitrary surfaces of revolution, but other less accurate methods of 
analysis have been developed for blades with arbitrary shape (refer- 
ence 7, for example) . The correlation equations to be developed are an 
extension of the correlation equations developed for straight blades on 
conic flow surfaces of revolution in reference 6. This extension enables 
the correlation equations to apply to arbitrary surfaces of revolution. 

In this report, as well as in reference 6, the correlation equations are 
developed to predict the distribution of the relative velocity components 
Qg and and the stream function \|r throughout the impeller passage. 

The correlation equations are developed from the fundamental conditions 
of continuity and absolute irrotational motion and fhctm certain observed 
results of the relaxation solutions. These correlation equations predict 
flow conditions in straight-bladed impellers with any tip speed, flow 
rate, area variation, blade spacing, and for any flow surface of revolu- 
tion. The correlation equations determine these flow conditions from 
various parameters that describe the design and operating conditions of 
the impeller and from tabulated values of any known "standard" solution 
(such as that given in reference 6) . 

In the correlation equations of this report, and of reference 6, the 
standard conditions (obtained from any known solution, reference 6, for 
example) and the nonstandard conditions (being investigated by the corre- 
lation equations) are evaluated at the same values of | and tj. Tbe 
effect of the shape of the surface of revolution on the fluid properties 
at a given point (|,tj) on the surface of revolution is determined by 
the effect of the shape on sin a, which term appears in the correlation 
equations, and by the effect of the shape on |, which quantity is 
related to sin a by equation (14a) . | 

Compressible flow. - In the previous development of correlation 
equations the local velocity was found to be, in effect, made up of the 
average value of at the given radius and superimposed velocity com- 

ponents caused by the rotation of the channel. The superimposed velocity 
components were sufficiently small not to be greatly affected by compress- 
ibility^ which was only considered to affect the average value of Q^. 

Because this effect on is easily incorporated, the correlation equa- 

tions in this report are developed for compressible as well as incompress- 
ible flow. 
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Assume -that the Q q component of the relative velocity is corre- 
lated hy 


Qe* =a(q 0 ) s ‘ ( 21 ) 

where the superscript prime indicates approximate values determined hy 
the correlation equations, where the subscript s indicates standard 
values (given in tables of reference 6, for example) , and where 

(AS) sin a 

A = [Mp (AS) sin aj g ^ 

For constant a equation (21) is identical to the corresponding corre- 
lation eqution of reference 6. (Note the different definition of AS 
in reference 6.) 

From equation (21) and using the identical procedure described in 
reference 6, the correlation equation for becomes 


Q m' “ ^Vav + A [Si " ^avjs + 

f(Q ) r 1 ^ 

A(2tj - 1) J ■ - j^(AS) sin a - (AS) g sin aj + (R - R s ) (Mp AS) s sin 

(23) 

where R is a known function of | and where from one- dimensional con- 
tinuity considerations (reference 6) 


Wav 



(24) 


In equation (24) the average density ratio P av /p Q is obtained from 

equation (9) by assuming that for straight blades the relative tangential 
velocity is zero and the meridional velocity Is equal to its average 
value so that 


Q 2 ~ Wav 2 


2466 



2466 


NACA TN 2654 


15 


Likewise the correlation equation for the stream function ijr becomes 



For a complete eddy flow (no through flow) * t and (Q M ) av are zero so 

that equation (25) becomes indeterminate. If, however, both sides of 
equation (25) are multiplied by 


*t - <Vav m 

which equality is obtained from equations (18) and (24), equation (25) 
can be evaluated directly for ijr' . 

jhcQmpressible flow. - For incompressible flow the stagnation speed 
of sbund can, as before, be eli min ated from consideration by expressing 
Q and ijr as ratios of 1^. Thus for incompressible flow the correla- 
tion equations become 




(21a) 


where 


a 


A* _ (A6) sin 

£(A0) sin a 


(22a) 






av 


% 


A* (2t) - 1) 


{i(4)s [ (A0) Sin a " (A0) s sin a s] + 

a ) 


(H - E s ) (A0) sin 


(23a) 
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and 



Repr es entat ive values of 



and 


♦i 


obtained by 


the relaxation solutions of this report are compared in figures 12 to 14 
with the corresponding values given by the correlation equations [(21a) , 
(23a) and (25a)] using the standard solution given in reference 6. 
Perfect agreement corresponds to the 45° straight line, and it is seen 
that the correlation equations are accurate. 


Effect of flow surface shape. - The correlation equation (23a) has 
been used to compute the velocity distribution on the blade surfaces for 
incompressible flow in impellers with straight blades, with constant 
flow area along the flow surface, and with various shapes of flow sur- 
faces of revolution. These flow surfaces are described by the shape of 
their generating lines (meridional streamlines) in the RZ-plane . In 
order to demonstrate the effect of flow surface shape, calculations have 
been carried out for generating lines that are described by the equation 



(26) 


The family of generating lines is given by various values for the 
parameter n, and as indicated in figure 15, the surfaqes in al l cases 
start at a radius R of 0.5 (impeller inlet) and end at a radius R of 
1.0 (impeller tip) . 

The resulting velocity distributions on the blade surfaces are 
plotted in figure 16. As the value of n decreases, the difference in 
velocity on the two blade surfaces toward the impeller tip decreases and, 
in the limit as n approaches zero the velocities become equal on both 
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to 

if* 

CD 

CD 


surfaces. At this limiting value of n equal zero the flow is axial 
(see generating line for n equal to zero in fig. 15) and the veloc- 


ities on both blade surfaces are equal to the average velocity 


(M 

\ MtMv’ 


which for these examples is constant because of the constant prescribed 
flow area. 


As the value of n increases the difference in velocity on the two 
blade surfaces toward the impeller tip increases and in the limit as n 
approaches infinity this difference approaches that for a radial flow 
impeller (see generating line for n equal to infinity in figure 15) . 

If the velocity distribution curves in figure 16 were extrapolated 
to lower values of R, the curves would cross because, 4 s indicated in 
figure 15, the flow surfaces that are more nearly radial near the 
impeller tip are more nearly axial near the impeller inlet, and vice 
'versa. 


SUMMARY OF REStffiES 

A general method of analysis is developed for Bteady, two- 
dimensional, compressible, nonviscous, irrotational flow between arbi- 
trary blade shapes on arbitrary flow surfaces of revolution in turbo- 
machines. Incompressible solutions are presented for four flow rates 
through each of two centrifugal impellers with the same hub-shroud con- 
tours but with different blade spacings. In addition, correlation equa- 
tions are developed whereby the velocity components and the stream func- 
tion distribution can be predicted lor compressible or incompressible 
flow in straight-bladed impellers only, with any tip speed, flow rate, 
area variation, blade spacing, and for any flow surface of revolution. 


Lewis Flight Propulsion Laboratory 

Rational Advisory Committee for Aeronautics 
Cleveland, Ohio, December 23, 1951 
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APPEKDIX 


SYMBOLS 


A 

A* 

arjr 


The following symbols are used in this report: 
parameter, equation (22) 
parameter, equation (22a) 
flow area at c han n el tip, equation (19a) 
stagnation speed of sound "upstream of turhamachine. 


tion (l) 


equa- 


g 

H 

h 

n 

P 

P 


gravitational acceleration 
height ratio, equation (10) 

spacing between center lines of channels farmed by streamlines 
adjacent to streamline that generates surface of revolution, 
figure 2 

blade-tip Mach number, equation (2) 
exponent, equation (26) 
pressure ratio, equation (8) 
static pressure 


Q relative velocity expressed as ratio of c Q 

meridional component of Q, equation (3) • 

QfijQfbQz components of Q in R, 0, and Z directions, respectively 

R,0,Z cylindrical coordinates, dimensionless (linear coordinates R 

and Z expressed as ratios of rrji) 


R" 


perfect gas constant 


rrp 


blade -tip radius 


T 


static temperature 


W 


total flow rate through channel (of height h) between two 
blades 
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a angle, figures 3 and. 4 

T dimensionless fluid circulation 

y ratio of specific heats 

T] transformed coordinate, equation (14b) 

Ad angular "blade spacing about axis of turbomachine 

X whirl ratio, equation (7) 

|i slip factor 

| transformed coordinate, equation (14a) 

p static weight density 

cp channel flow coefficient, equation (19) 

ijr stream function, equation (13) 

as angular velocity of blade row 

Subscripts: 

av average 

i incompressible 

min minimum 

0 stagnation condition -upstream of turbomachine 

s standard solution 

T blade tip 

t trailing face of blade 

« 

U far upstream of turbomachine 

Superscripts: 

1 estimated by correlation equations 
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TABLE I - CHARACTERISTICS OF MEAN FLOW SURFACE OF 
REVOLUTION USED I N NUMERICAL EXAMPLES 


i • Ad 

R 

z 

sin a 

H 

(Qg/MrjOu 

(a) 

-1.792 

0.4915 

-0.3394 

0.0082 

3.659 

1.001 

-1.664 

.4920 

- .2767 

,0140 

3.652 

1.002 

-1.536 

.4930 

- .2135 

.0218 

3.636 

1.004 

-1.408 

.4950 

- .1502 

.0333 

3.608 

1.008 

-1.280 

.4978 

- .0870 

.0555 

3.565 

1.015 

-1.152 

.5021 

- .0232 

.0900 

3.493 

1.027 

-1.024 

.5096 

.0411 

.1422 

3.364 

1.050 

- .896 

.5212 

.1058 

.2222 

3.160 

1.093 

- .768 

.5408 

.1710 

.3550 

2.858 

1.165 

- .640 

.5722 

.2348 

.5372 

2.544 

1.237 

- .512 

.6205 

.2933 

.7312 

2.236 

1.298 

. - .384 

.6890 

.3408 

.8888 

1.901 

1.374 

- .256 

.7770 

.3737 

.9678 

1.542 

1.502 

- .128 

.8803 

.3931 

.9913 

1.239 

1.651 

0 

1.0000 

.4052 

.9968 

1.0000 

1.800 

.128 

1.1368 

.4120 

.9996 

.8700 

1.820 

.256 

1.2919 

.4127 

1.0000 

.8664 

1.608 


a Values of ( Qg/Mj ) U f° r flew rates used in 
numerical examples of tils report are given 
in the text under Operating conditions. 
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Blade tip 


Shroud - 


Stream function 
ratio 

1.0 

.9 


Mean stre amlin e used 
to generate surface 
of revolution for 
solutions of this y 
report -i / 



Axis of turhomachine 


Figure 2. - Streamlines in meridional plane for axial-symmetry 
solution of channel described in figure 7. Streamline 
designation Indicates percentage of flow through channel 
between streamline and hub. Incompressible flow; Q^/Mp 
equal to 0.3429. far upstream (reference 4). 


Axis of 
turi) etna chine 





Figure 3. - Cylindrical coordinates and velocity components relative to rotating "blade 
row. All quantities are dimensionless. Linear coordinates are measured in uni ts of 
blade tip radius ; velocity components are measured in units of stagnation speed of 
sound upstream of blade rov. 
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Figure 4. - Coordinates, velocities, and so forth 
at point on streamline in meridional plane (axial- 
symmetry solution) . 



Front view Side view 


Figure 5. - Fluid particle on developed surface of cone at 
point that-ls' tangent tro flow surface of revolution. 
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co 

co 

cvi 



Axial distance, Z 



Figure 7. - Hub-shroud dimensions of impeller for numerical 
examples . Vaneless diffuser; straight imp eller blades 
extended Indefinitely far upstream parallel -with- axis -of 
Impeller. 
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(b) Meridional velocity component. 


% 

*P 


Figure 11. - Conoluded. Lines of constant relative 
velocity cocaponenta. AQ, 14.67*; eddy flow only. 
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Figure 14. - Comparison between relaxation and estimated values of ifa. 
Correlation equation (25a). 
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.5 .6 .7 .8 .9 1.0 


Axial distance ratio, Z 

Figure 15. - Generating lines in axial-radial plane that are 
rotated about axis of impeller to obtain flow surfaces of 
revolution. 











